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We report on the synthesis, STM imaging and theoretical studies of the structure, electronic 
structure and transport properties of linear chains of styrene and methyl-styrene molecules and 
their heterojunctions on hydrogen-terminated dimerized silicon (001) surfaces. The theory presented 
here accounts for the essential features of the experimental STM data including the nature of the 
corrugation observed along the molecular chains and the pronounced changes in the contrast between 
the styrene and methyl-styrene parts of the molecular chains that are observed as the applied bias is 
varied. The observed evolution with applied bias of the STM profiles near the ends of the molecular 
chains is also explained. Calculations are also presented of electron transport along styrene linear 
chains adsorbed on the silicon surface at energies in the vicinity of the molecular HOMO and LUMO 
levels. For short styrene chains this lateral transport is found to be due primarily to direct electron 
transmission from molecule to molecule rather than through the silicon substrate, especially in the 
molecular LUMO band. Differences between the calculated position-dependences of the STM current 
around a junction of styrene and methyl-styrene molecular chains under positive and negative tip 
bias are related to the nature of lateral electron transmission along the molecular chains and to the 
formation in the LUMO band of an electronic state localized around the heterojunction. 
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I. INTRODUCTION 



During the last few years there has been growing in- 
terest in molecular electronics, stimulated largely by the 
experimental realization of molecular wires^*2^ systems 
in which a single organic molecule or a few molecules 
carry an electric current between a pair of metal contacts. 
In some cases such systems exhibit switching behavior 
and/or negative differential resistances^ phenomena 
that may be exploited in future molecular electronic de- 
vices. Hybrid molecular/semiconductor nano-electronic 
devices are another intriguing possibility and fundamen- 
tal research that may ultimately lead to their creation is 
also being pursued at the present timeASaLiSiiiiiLiiii^ 
Recently, techniques have been developed that make 
it possible to grow an oriented linear chain of styrene 
molecules on a hydrogen-terminated silicon substrate be- 
ginning at a predefined point on the substrate The 
electronic structure and STM images of these homo- 
molecular chains of styrene have also been modelled 
theoretically^^ However, no studies (experimental or 
theoretical) of chains of styrene molecules with chemical 
substituents attached have been reported to date. In the 
work reported in the present article we extend the exper- 
iments and theory to a related system, a heteromolecular 
chain of molecules on a silicon substrate where the chem- 
ical identity of the molecular species changes abruptly 
across a heterojunction from styrene to methyl-styrene. 
As the chemical composition across the molecular chain 



varies solely by the presence (or absence) of a substituent 
(in this case a methyl group) on the aromatic ring, the tt- 
7t* stacking that is an essential prerequisite for wire-like 
conduction along the chain remains undisrupted. 

Some intriguing questions regarding the self-assembled 
lines of molecules on silicon are: To what extent do 
the molecules that form the line communicate with each 
other electronically? Is electrical conduction along the 
line possible, and if so is such conduction dominated by 
transport through the molecules themselves or through 
the silicon substrate? Calculations of lateral trans- 
port along the molecular chain can in principle answer 
these questions. However, to date such calculations 
have not been reported for molecular chains on sili- 
con, although molecular band formation in the styrene 
chains due to orbital overlap between molecules has been 
studied)2iiS as has conduction through stacks of aromatic 
biphenyldithiol molecules not attached to a substrate^ 
The role of molecule-molecule interactions in vertical con- 
duction through molecules between metal contacts has 
also been investigated theoreticallyiiSiiLiSiiSiSSi 

To date it has not been feasible experimentally to 
connect a pair of probes to the two ends of a chain of 
molecules on silicon so as to measure the lateral trans- 
port directly: STM studies of the molecular chains mea- 
sure vertical transport between the silicon substrate and 
a single metal tip via the adsorbed moleculesi&iS How- 
ever, vertical transport measurements may yield indi- 
rect experimental information relevant to lateral trans- 
port, particularly if the molecular chain is inhomoge- 
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neous. A heterojunction between chains of similar but 
distinct molecules such as styrene and methyl-styrene on 
silicon or the end of a line of molecules may be a suitable 
inhomogeneity for such studies, but theoretical model- 
ing of such systems is clearly needed in order to extract 
meaningful information of this kind from experimental 
data. 

A realistic model is presented here of the structure, 
electronic structure and electronic transport in chains 
of styrene and methyl-styrene molecules on hydrogen- 
terminated dimerized silicon (001) surfaces. The model 
accounts for the main properties of heterojunctions of 
styrene and methyl-styrene molecular chains on silicon 
that have been observed in our STM experiments, as 
well as the experimentally observed phenomena associ- 
ated with the ends the molecular lines. The first calcu- 
lations are presented of lateral transport though short 
styrene chains adsorbed on the silicon surface. It is 
found that lateral electron transport along an 8-molecule 
styrene chain at energies corresponding to the molecu- 
lar HOMO and LUMO bands is primarily due to direct 
electron transmission from molecule to molecule rather 
than transmission via the silicon substrate. This dom- 
inance of molecule-to-molecule transmission over trans- 
mission through the substrate is stronger for the LUMO 
band than the HOMO band. This difference is related to 
pronounced differences between the calculated position- 
dependences of the STM (vertical transport) current 
around the junction of styrene and methyl-styrene molec- 
ular chains under positive and negative tip bias. 

This article is organized as follows: In Section ITU we 
present experimental results detailing the growth and 
STM imaging characteristics of heteromolecular chains of 
styrene and methyl-styrene on silicon. In Sections IIIII and 
IIVI we describe our modeling of the structure and elec- 
tronic structure of the molecular chains. The approach 
used in our transport calculations is summarized in Sec- 
tion [V] Our theoretical results for vertical transport in 
the STM geometry in the vicinity of a heterojunction 
of molecular lines are presented and compared with ex- 
periment in Section lvTl Our theoretical results regarding 
lateral transport along the molecular chains and how this 
may influence STM images are presented in Section [VIII 
Our experimental and theoretical findings regarding the 
structure and STM images of the ends of the molecular 
chains are discussed in Section lYlIII Finally, our conclu- 
sions are summarized in Section llXI 



II. EXPERIMENT 

The experiments were performed on a hydrogen termi- 
nated Si(100)-2xl surface. Samples were cleaved from 
an arsenic-doped (0.005 f2-cm) Si (100) oriented wafer, 
mounted in molybdenum sample holders, and load- locked 
into an ultrahigh vacuum (UHV) chamber (background 
system pressure < lxlO -10 Torr). Samples were de- 
gassed for 8 hours at 700°C. Removal of the surface oxide, 



and sublimation/replanarisation of the surface layers was 
achieved by repeated thermal cycling to 1250°C. Sam- 
ple temperatures were monitored using a calibrated in- 
frared pyrometer. During the annealing procedure, sam- 
ple heating was discontinued if at any point the mon- 
itored system pressure exceeded 4xl0~ 10 Torr. This 
sample preparation method produces 100 oriented silicon 
surfaces with well ordered silicon terraces (2x1 surface 
reconstruction), and a surface defect concentration below 
5%. 

The vacuum preparation chamber was then filled with 
hydrogen gas (lxlO -6 Torr), and a hot tungsten filament 
(T ~ 1600°C) situated 10 cm from the sample was used to 
provide a flux of reactive atomic hydrogen to the surface. 
The sample was heated to 330° C and exposed to the 
atomic hydrogen flux for 13 minutes. These conditions 
allowed the formation of a quasi-saturated monolayer of 
silicon monohydride with a 2x 1 surface reconstruction. 21 

One-dimensional molecular organic heterowires were 
grown along dimer rows of the as-prepared H:Si (100) 
surfaces using a vacuum phase reaction studied earlier by 
Lopinksi and co-workers £ In the case of the styrene reac- 
tion with the H:Si surface, the double carbon bond in the 
vinyl group at the base of styrene molecule reacts with a 
surface silicon radical (corresponding to a hydrogen va- 
cancy on the H:Si surface) to produce a silicon-carbon 
bond and a carbon radical (centred on the carbon atom 
bound to the aromatic ring). The carbon centred rad- 
ical then abstracts a hydrogen atom from an adjacent 
dimer on the same row (and from the same side of the 
dimer - never from the other (diagonal) atom), produc- 
ing a newly reactive site. Successive reactions lead to the 
well ordered linear structures reported in Ref. |8j. Sev- 
eral substituted styrene compounds have been studied in 
our group and have demonstrated this self-directed line 
growth mechanism on H:silicon^ Here two such com- 
pounds are used, styrene and 4-methyl-styrene, to fabri- 
cate one-dimensional organic heterowires. 

On the H:Si(100) surface, residual silicon radicals typ- 
ically exist at the level of a few percent and are used in 
this work as initiation sites for the self-directed growth 
of 4-mcthyl-styrene/styrene heterowires. The molecu- 
lar heterowires were fabricated by sequentially introduc- 
ing 4-methyl-styrene and styrene into the STM cham- 
ber through a variable leak valve. Dissolved atmospheric 
gases in the styrene and 4-methyl-styrene were expelled 
from the liquid phase samples by repeated freeze-pump- 
thaw cycles. Prior to dosing, the STM tip was retracted 
~ 1 micrometer from the H: silicon surface. During dos- 
ing, the flow of vapour phase reactant was adjusted to 
bring the chamber pressure to ~ 4xl0~ 7 Torr, and held 
for several tens of seconds (integrated doses given below) . 
After dosing (valve shut) , the chamber pressure returned 
below lxl0~ 10 Torr. 

Images of the H:Si(100) surfaces were collected in UHV 
at room temperature using an Omicron STM1. Electro- 
chemically etched tungsten tips were cleaned in-situ by 
electron beam bombardment and field emission prior to 
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STM imaging. Constant current topographic STM im- 
ages were acquired using a constant tunnel current of 
60 pA. Quoted dimensions along the lateral and vertical 
directions have been scaled to agree with known values 
for the silicon dimer spacing, and silicon terrace height, 
respectively. 4-methyl-styrene/styrene heterostructures 
were grown on several H-silicon surfaces, and studied 
with different STM tips. 

Figs. and show the growth and bias-dependent 
filled-state imaging of a methyl-styrene/styrene het- 
erostructure. Fig^a), acquired using a sample volt- 
age (V s ) of -3.0 V, shows the formation of the first line 
segment following a 8 Langmuir dose of methyl-styrene 
(1 Langmuir = 10~ 6 Torr-seconds). Approximately 8 
molecules of methyl-styrene have reacted at the indicated 
location. The bright feature at the end of the line seg- 
ment (indicated by the white arrow) corresponds to the 
terminal silicon radical (i.e. a bare silicon atom which 
has transferred its hydrogen atom to the adjacent methyl- 
styrene molecule). Its presence indicates that the reac- 
tion between the methyl-styrene and the silicon surface 
at this particular site has not reached saturation. 

Figdb) shows the same region of the crystal following 
a 78 Langmuir exposure of styrene. Starting at the silicon 
radical in (a) , approximately 8 molecules of styrene have 
subsequently reacted to form a methyl-styrene/styrene 
heterostructure. The location of the heterointerface is 
indicated by the white wedge. The silicon radical present 
in (a) is absent from the end of the styrene line segment 
in panel (b), and in panels (c) and (d) acquired at lower 
magnitude bias. This suggests that linegrowth has termi- 
nated at an unresolved surface defect (such as a dihydride 
site), or has been capped by a H-atom deposited by the 
STMtip^ 

Panels (b) - (d) show the topographic envelope of the 
heterostructure at varying sample bias. At V s = -3.0 
V [panel (b)], the constant current contours follow a 
relatively smooth envelope modulated by the physical 
height of the methyl-styrene and styrene line segments 
above the FLsilicon surface. At this value of the bias the 
aromatic it states in the styrene and methyl-styrene lie 
above the tip Fermi-level. The tunnel current therefore 
results from electrons which tunnel from underlying oc- 
cupied states within both the silicon and surface bound 
molecules into empty states within the STM tip. The ex- 
tent to which molecular orbital states emptied by tunnel- 
ing are replenished directly by underlying silicon states 
(i.e. vertical transport across molecules) or indirectly by 
nearest neighbour molecules (reflecting a degree of lat- 
eral transport across the molecular chain) is considered 
in the following theory sections. 

In panel (c) [V s = -2.4 V], the aromatic it states move 
below the tip Fermi-level, and the imaging current is sup- 
plied by silicon states which tunnel across the surface 
bound molecules. In this regime, the methyl-styrene and 
styrene image with similar height. Also, molecular scale 
corrugation appears within the topographic envelope. In 
panel (d) [V s = -1.8 V], the molecular corrugation shows 



different contrast between the methyl-styrene and styrene 
line segments. 

To better resolve differences in the imaging character- 
istics between both segments, 0.2 nm wide topographic 
cross-sections were extracted along two parallel axes 
within the heterostructure envelope. Images were reg- 
istered using a least squares fitting algorithm to allow 
extraction of cross-sections from nominally identical lo- 
cations between images acquired at different bias. 

Cross-sections taken along the direction of the blue 
arrow (Fig^ panel (d) inset) appear in panel (a) of Fig[21 
and are centred above the chemical attachment points 
between the heterostructure and the underlying dimer 
row. Cross-sections taken along the direction of the red 
arrow, appear in panel (b) of Fig[2 and are centred to 
the right of the trench separating the underlying dimer 
rows. Cross-sectional data extracted from images (not 
shown) acquired at sample voltages of -2.0 V, -2.2 V, - 
2.6 V, and -2.8 V have been added to the data presented 
in FigHJ 

In panels (a) and (b) of FigEl the height envelope cor- 
responding to the methyl-styrene line segment is situated 
between 2.5 nm and 6.0 nm along the distance axis. The 
styrene line segment is situated between 6.0 nm and 9.8 
nm. The peaks centred at 1.5 nm (in both panels) and 
at 10.3 nm (in panel (a) only) correspond to the isolated 
(and unidentified) structures present in the upper and 
lower left hand corners of the images in Fig^ 

In each panel the position of the corrugation maximum 
associated with the interfacial methyl-styrene molecule is 
indicated by a black arrow. The position of the chemical 
interface in the heterowire is made evident by the dif- 
ferential height contrast observed between the two line 
segments at high and low bias. Superposition of the het- 
erowire envelope with registered cross-sections extracted 
from Fig^ (a) (not shown) agree with this determina- 
tion. In this case, the topographic maximum (i.e. the sil- 
icon radical) at the end of the methyl-styrene segment in 
Figffl(a) lies between the as-identified interfacial methyl- 
styrene and styrene corrugation maxima in the com- 
pletely grown heterowire. Other methyl-styrene/styrene 
heterowires studied displayed similar agreement between 
the apparent and expected locations of the heterojunc- 
tion. 

The main features displayed in Fig^ are visible in 
the cross-sectional data presented in panels (a) and (b) 
of Fig [3 As the magnitude of the imaging bias is in- 
creased, the apparent height of the methyl-styrene and 
styrene line segments increases and begins to saturate 
(with methyl-styrene imaging above styrene). At low 
magnitude bias (where the contribution of molecular it 
states to the tunnel current is suppressed) , the molecular 
corrugation is more apparent. 

While in panel (b) the height of the methyl-styrene 
line segment is comparable to (at low magnitude bias) 
or greater than (at high bias) that of styrene, a dif- 
ferent behaviour is observed along the chemical attach- 
ment points in (a). At low magnitude bias, \V\ < 2.6 
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V, the methyl-styrene molecules image lower than the 
styrene molecules. Given that methyl-styrene extends 
farther beyond the H:silicon surface than styrene, and 
the similarity in molecular conformation anticipated for 
the molecules in each segment, this observation seems 
unexpected. 

At the heterojunction, the imaging contrast between 
the interfacial molecules (styrene and methyl-styrene) ap- 
pears to a first approximation comparable to the contrast 
resolved for these same molecules within their respec- 
tive line segments. A significant departure from this be- 
haviour is observed at elevated magnitude bias (\V\ ^ 2.6 
V), where the interfacial methyl-styrene appears lower 
than the molecules located within the line segment. The 
molecules at either end of heterowire also image with de- 
creased height. Finally, a slight upward bowing (roughly 
centred at the heterointerface) appears superimposed on 
the bulk of the cross-sectional curves. 

In the absence of electronic interaction between 
molecules and geometrical differences between them, ev- 
ery molecule along a homomolecular wire (or a homo- 
molecular line segment in the case of a heterowire) would 
image with identical properties. This is not observed. 
The implication is that the observed height profiles along 
the heterowire reflect some combination of: i) conforma- 
tional changes in the molecules along the chain due to 
intermolecular forces, ii) dispersion interactions leading 
to differential broadening and shifting of molecular elec- 
tronic levels along the chain, and/or iii) tunneling current 
collected by the STM tip above a given molecule receiv- 
ing contributions from neighbouring molecules along the 
chain reflecting some degree of lateral current transport 
across the stacked aromatic rings. 

To determine the origin of the observed contrast in 
molecular corrugation, and the topographic height varia- 
tions resolved along the heterowire structure, a theoreti- 
cal model for the heterowire is developed and its current 
transport properties are considered in sections IV to VIII. 



III. MODELING THE STRUCTURE 

The spacing between the tops of the molecules seen in 
our experimental STM images of styrene/methyl-styrene 
chains on the (001) silicon surface varies along the length 
of the chain and is larger near the ends of the chain 
than near its center. Because of the variable mismatch 
between the molecular spacing and the underlying sili- 
con lattice, it is not appropriate to model these molecu- 
lar chains theoretically as structures periodic along the 
length of the chain with a small unit cell. They are there- 
fore modelled here as finite clusters such as that depicted 
in Fig|3]which shows a molecular chain with a heterojunc- 
tion, consisting of four styrene and four methyl-styrene 
molecules. The molecules bond to a dimerized Si (001) 
surface that is represented by two rows each of 10 sili- 
con dimers together with four underlying atomic layers 
of silicon atoms. All of the dangling bonds are passivated 



with hydrogen. The lowest carbon atom of each molecule 
bonds to a surface Si atom and the molecules lean over 
the "trench" between the two rows of silicon dimersiS^ 

A complete ab initio relaxation of this cluster using 
Kohn-Sham density functional theory would be imprac- 
tical, so the following hybrid approach was adopted to 
calculate its approximate structure: An ab initio den- 
sity functional relaxation was carried out^l of a smaller 
cluster consisting of one styrene and one methyl-styrene 
molecule over two rows each of two silicon dimers, with 
two additional underlying layers of silicon, all of the dan- 
gling bonds passivated with hydrogen. This small clus- 
ter shown in Figgis similar to the immediate vicinity of 
the junction between styrene and methyl-styrene chains 
in Fig|3] The relaxation of this small cluster was carried 
out keeping the atoms of the two lower silicon layers fixed 
at positions corresponding to a bulk silicon crystal lat- 
tice but allowing the positions of all of the atoms of the 
molecules, the silicon dimer atoms and all of the hydro- 
gen atoms to relax freely. The main difference between 
the relaxed geometries of the two molecules in FigQ] and 
that of an isolated molecule bound to the same silicon site 
is that the two molecules tilt somewhat away from each 
other. This tilt occurs mainly through changes in the 
sequence of dihedral angles defined (in the Z-matrix de- 
scription of the system) by the orientations of the bonds 
connecting Si atom 1 to C atom 2, C atom 2 to C atom 3, 
and so on along the carbon atom chain through C atom 
6 of the styrene molecule (see FigQJ, and of the cor- 
responding bonds of the methyl-styrene molecule. Thus 
the relaxed geometry of the cluster in Fig|3]was generated 
by assigning initially to all of the atoms of the molecules 
positions derived from the relaxed geometries of the re- 
spective atoms of the smaller cluster in FigQJ Then this 
initial structure was relaxed so as to minimize the en- 
ergy of the cluster using the AMI empirical model^ but 
allowing only the above dihedral angles, the dihedral an- 
gles defining the orientations of the C-H bonds involv- 
ing carbon atoms 2 and 3 and corresponding C-H bonds 
on other molecules, and the dihedral angles defining the 
orientation of the triplets of H atoms belonging to the 
methyl groups to vary. During the AMI relaxation a fur- 
ther constraint was also imposed, namely, that the tilts 
of the styrene and methyl-styrene molecules at the ends 
of the 8 molecule chain be such that the average lateral 
spacing between the top carbon atoms of the molecules 
approximate the experimentally observed average spac- 
ing of the molecules near the junction of the styrene and 
methyl-styrene chains^ The final model geometry of the 
8 molecule chain is thus representative of the region of a 
much longer molecular chain around the heterojunction. 
A similar procedure was used to obtain the relaxed ge- 
ometry of the chain of 8 styrene molecules used in the 
theoretical studies of lateral electron transport along the 
molecular chain that are discussed in Section lYlII How- 
ever, the structures used in our theoretical studies of the 
ends of the molecular chains were obtained somewhat 
differently as is described in Section IVIIII Finally, it 
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was found that the energy of the structure with methyl- 
styrene molecules varies by much less than kT (where 
T = room temperature) when the methyl groups are ro- 
tated about the bonds connecting them to the benzene 
rings of the molecules. Since the experiments on these 
systems were carried out at room temperature, the cal- 
culated STM tip currents presented below are averages 
over orientations of the methyl groups. 

The atomic structure of the tungsten STM tip in our 
experiments is unknown and the tip was modelled arbi- 
trarily as a (OOl)-oriented (relative to the (001) Si sur- 
face) clean bcc tungsten tip terminating in a single tung- 
sten atom. A total of 15 tungsten atoms arranged in 
layers of 1,4,5,4 and 1 atoms along the (001) direction 
were included in the model of the tip. These atoms were 
placed at their nominal positions in a bulk bcc tungsten 
lattice except for the terminal atom whose position was 
relaxed along the (001) direction (towards the center of 
the tungsten cluster) so as to minimize the Kohn-Sham 
energy of the 15 atom tungsten cluster^ The effect of an 
alternate tip on the calculated STM currents, obtained 
by replacing the apex atom with a single silicon atom 
whose position was relaxed appropriately, was also inves- 
tigated. 



IV. MODELING THE ELECTRONIC 
STRUCTURE 

Most theoretical work on electronic transport in molec- 
ular wires with metal contacts has been based on 
semi-empirical tight-binding models or ab initio Kohn- 
Sham density functional calculations of the electronic 
structurei 3 i 31 i 32 i 33 i 34 i 35 i 36 However, the atomic structures 
of the molecule-metal junctions are not known. Thus, 
with a few exceptions, 36,37 comparison of theoretical 
predictions with experiments on metal/molecule/metal 
wires is subject to large uncertainties. By contrast, 
the bonding geometries of many small organic molecules 
to silicon substrates are known^ as are the specific 
atomic silicon sites involved in the carbon-silicon bonds 
anchoring styrene and methyl-styrene molecular chains 
to the silicon .SaLiiLii Thus comparison between theory 
and experiment is subject to much less uncertainty for 
these molecules on silicon than for metal/molecule/metal 
wires. This allows us to assess the validity of different 
models of the electronic structure of the molecular chains 
on silicon by comparison with experimental data. 

Previous theoretical work modeling the current-voltage 
characteristics of single styrene molecules on silicon found 
ab initio Kohn-Sham density functional calculations un- 
suitable for treating the electronic structure of the silicon 
substrate^ 3 - The present work on styrene and methyl- 
styrene chains on silicon reached a similar conclusion: 
Our density functional calculations that employed local- 
ized orbital bases appropriate for quantum chemistry^, 
in common with previous density functional calculations 
for styrene chains on silicon that employed plane wave 



basesfi - yielded molecular HOMO levels aligned with the 
highest occupied silicon valence states. This level align- 
ment is unrealistic; experimentally the molecular HOMO 
energy is considerably lower. Furthermore, if it were cor- 
rect then STM images of the styrene molecules probing 
occupied states at low bias would show a current node at 
the center of each molecule since the molecular HOMO 
is a 7r-like orbital on the benzene ring. However, this 
is not seen experimentally. Finally, density functional 
calculations in the local density approximation seriously 
underestimate the band gap of bulk silicon. For these 
reasons 38 a different approach based on extended Hiickel 
theory will be adopted here in modeling the electronic 
structure of the molecules, the silicon substrate and the 
tungsten STM tip. Our model corrects the above de- 
ficiencies of density functional theory and yields current 
maxima over the centers of the molecules for low negative 
applied substrate bias, consistent with the nature of the 
corrugation along the molecular chain that is observed 
experimentally. 

Extended Hiickel theory is a semi-empirical tight bind- 
ing scheme from quantum chemistry that provides an ap- 
proximate description of the electronic structure of many 
molecules. It describes molecular systems in terms of 
a small set of Slater-type atomic orbitals their 
overlaps Sij = (<^>i|</>j) and a Hamiltonian matrix Hij = 
(<t>i\H\(j)j). The diagonal Hamiltonian elements Ha = ti 
are chosen to be the atomic orbital ionization energies. 
In the Wolfsberg-Helmholz form of the model, the non- 
diagonal elements are approximated by 

H ij = KS ij (e i + e j )/2 (1) 

where K is a phenomenological parameter usually chosen 
to be 1.75 for consistency with experimental molecular 
electronic structure data. The extended Hiickel model 
has been used successfully to explain the experimen- 
tal current-voltage characteristics of a variety molecular 
wires connecting metal electrodes^S* 3 ^ However, with- 
out modification it is unsatisfactory for crystalline silicon 
as it predicts a band gap that is direct and roughly three 
times larger than that found experimentally for this ma- 
terial. In the present work these deficiencies of the ex- 
tended Hiickel model have been corrected by replacing 
equation Q for the Hamiltonian matrix elements be- 
tween the atomic orbitals of silicon by 

H tj = K i0 - (ei + 6^/2 (2) 

where the Kij are fitting parameters*! chosen so that 
the modified extended Hiickel model obtained in this way 
yields the correct silicon band structure. The extended 
Hiickel Hamiltonian matrix elements for tungsten were 
adjusted similarly so as to match the band structure of 
bcc tungstenjai however in the case of tungsten some 
modification of the extended Hiickel model's orbital en- 
ergies was also found to be necessary to achieve a good 
fit. Finally a realistic alignment between the molecu- 
lar HOMO level and the Fermi levels of the silicon and 
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tungsten clusters was achieved by shifting the silicon and 
tungsten orbital energies by appropriate amounts Egi 
and Ew respectively^ Because of the non-orthogonality 
of the extended Huckel basis states, such shifts of the di- 
agonal matrix elements of the Hamiltonian also required 
corresponding shifts AHij of the non-diagonal elements 
Hij , which were taken to be 

(3) 

where Ek — Esi or Ew if {| </>&}} is a silicon or tungsten 
orbital, respectively. 

FigEJshows the key features of the electronic structure 
of the eight-molecule styrene/methyl-styrene heterojunc- 
tion on the hydrogen passivated silicon cluster whose ge- 
ometry is depicted in Fig|3 as described by the above 
electronic structure model: The energy eigenvalues and 
eigenstates of the model Hamiltonian were calculated and 
a Mulliken analysis was carried out to determine the sil- 
icon, carbon and hydrogen content of each eigenstate. 
The carbon content was further resolved according to 
whether the carbon atom belongs to a styrene or methyl- 
stryrene molecule. The results are shown as histograms 
in FigOa (Si), b (C/styrene), c (C/methyl-styrene) and d 
(H), summed for the sake of clarity over the eigenstates in 
bins of energy of width 0.1 eV. The vertical line labelled 
"Si HOMO" marks the energy of the highest occupied 
state of the cluster which has mainly silicon content. 
The carbon band below about -12.3 eV corresponds to 
molecular HOMO orbitals while the carbon band around 
-8.2 eV corresponds to the molecular LUMO. Because the 
HOMO (LUMO) of ethylbenzene is close in energy and 
generally similar to the HOMO (LUMO) of methylethyl- 
benzene, the molecular HOMO (LUMO) bands due to 
the styrene and methyl-styrene on silicon (Fig[5J} and c) 
are similar to each other although the styrene band edges 
are slightly lower (by ~ 0.1-0.2 eV) in energy than those 
for methyl-styrene. 

Note also the weak carbon features between the molec- 
ular HOMO band edge at -12.3 eV and the (mainly sili- 
con) cluster HOMO at -11.3 eV. These are almost iden- 
tical for the styrene and methyl-styrene molecules and 
are due to electron tunneling through the molecules from 
occupied states of the silicon substrate that lie above 
the molecular HOMO band. This tunneling is respon- 
sible for the STM current for low negative sample bias 
at which the experimentally observed STM images of 
the molecules are similar to those calculated using the 
present model (Fig|Bj plot L in Section IVlf) but differ 
markedly from the those predicted by density functional 
calculations as was discussed above. 

Since the styrene and methyl-styrene molecules are 
coupled more strongly to the silicon substrate than to 
the STM tip, it is assumed for simplicity in the present 
work that when a bias voltage V is applied between the 
STM tip and sample, the entire voltage drop occurs be- 
tween the molecules and the STM tip4£ The correspond- 
ing shift eV of the energy levels of the tip relative to the 
sample is then included in Ew, and thus also contributes 



to the off-diagonal matrix elements of the Hamiltonian 
through equation 

Finally, it should be emphasized that because the 
present model includes only a few layers of silicon atoms, 
the electronic structure of the silicon that it describes is 
characteristic of the immediate vicinity of the silicon sur- 
face only and does not include band bending effects that 
occur over larger length scales within the silicon. Fur- 
thermore the bias voltages considered in the model are 
those between the region slightly below the silicon sur- 
face and the STM tip, whereas bias voltages measured 
experimentally are those between the STM tip and the 
deep silicon interior that is separated from the surface 
by regions where significant bias-dependent band bend- 
ing occurs. Thus from the stand point of theory the ap- 
plied bias will be characterized according to the location 
of the Fermi level of the tungsten tip (or, more precisely, 
its electrochemical potential) relative to the features of 
the carbon and silicon partial densities of states that are 
shown in Fig[3 rather than the numerical value of exper- 
imental bias voltage. 

V. MODELING THE ELECTRONIC 
TRANSPORT 

In the transport calculations reported here the elec- 
tronic structures of the silicon substrate, molecules and 
tungsten tip are all treated together as a single system de- 
scribed by a single Hamiltonian. Thus, unlike in STM im- 
age calculations based on Bardeen tunneling theory, here 
the Hamiltonian matrix elements connecting orbitals of 
the molecules and those of the tip are not assumed to 
be small. Because of this the STM current can be calcu- 
lated for small separations between the tip and molecules 
as well as large. 

The calculations of the current are based on Landauer 
theoryM which relates the electric current / through a 
nanostructure under an applied bias voltage V to the 
multichannel electron transmission probability T(E, V) 
through the nanostructure at energy E asa^ 

2p I" 00 

W = Xi dET(E,V)(f{E,n a )-f{E,n d )) (4) 

where f(E, fii) = l/(exp[(i? — fj,i)/kT] + 1) and is the 
electrochemical potential of the source [i = s) or drain 
(i = d) electrode^ 

In the present work the source and drain electrodes 
are modelled as arrays of ideal leads: For calculations 
of vertical transport in the STM geometry, the electrode 
connected to the silicon substrate is modeled as an array 
of 152 ideal single-channel leads, one such lead coupled 
to the Is orbital of each of the hydrogen atoms that pas- 
sivate the silicon dangling bonds of the three lower layers 
of silicon atoms of the silicon cluster that represents the 
silicon substrate in Fig|3] The electrode connected to 
the STM tip was modelled as an array of 81 ideal single- 
channel tight-binding leads, one such lead coupled to each 
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of the 9 (s,p, d) valence orbitals of each atom of the model 
tungsten tip that is not the tip atom itself or one of its 5 
closest neighbors. As well as mimicking macroscopic elec- 
trodes by supplying an ample electron flux to the system, 
this large number of ideal source and drain leads has a 
similar effect to phase-randomizing Biittiker probes^ in 
minimizing the influence of dimensional resonances due 
to the finite sizes of the tungsten and silicon clusters em- 
ployed in the model. Each ideal lead i was modelled as a 
semi-infinite tight-binding chain with a single orbital per 
site, a site energy ct^ and nearest neighbor hopping ma- 
trix element f3. on was chosen to be equal to the energy 
(including the shift, if any, due applied bias described in 
Section llVjl of the hydrogen or tungsten orbital to which 
lead i was coupled. The hopping matrix elements were 
all taken to have a magnitude of 5 eV, sufficiently large 
that the eigenmodes of all of the ideal leads in the energy 
window between fj, s and fid in eq. be propagating;^, 
The coupling matrix elements Wj between the ideal leads 
and their respective hydrogen and tungsten orbitals were 
also set equal to (3. 

To calculate T(E,V) and hence evaluate eq.JU, the 
transformation to a different Hilbert space described in 
Refs. l48l and 1491 was made, mapping the non-orthogonal 
basis of atomic orbitals discussed in Section IIVI to an 
orthogonal basis. The Lippmann-Schwinger equation 



|*i> = \§o,i) + G a (E)W\*i) 



(5) 



describing electron scattering between the source and 
drain leads via the tungsten tip, molecules and hydrogen 
terminated silicon cluster was solved numerically for 
in the alternate Hilbert space42^ In eq.JSJ) G (E) is the 
Green's function for the decoupled system (i.e., with the 
coupling between the ideal leads and the hydrogen and 
tungsten orbitals switched off), is the eigenstate of 

the decoupled ideal source lead i with energy E and 
is the corresponding scattering eigenstate of the complete 
system with the coupling W between the ideal leads and 
the hydrogen and tungsten orbitals switched on. The 
scattering amplitudes tji from the ideal source lead i to 
drain lead j at energy E were extracted from the scat- 
tering eigenstates and the transmission probability 
that enters eq. was then calculated from 



(6) 



where V{ and Vj are the electron velocities in ideal leads 
i and j respectively at energy E. 

Lateral transport through chains of styrene molecules 
on the silicon substrate was modelled similarly. How- 
ever, in this case no STM tip was included in the model 
and no ideal leads were coupled to the silicon cluster or 
to the hydrogen atoms passivating it. Instead, 18 ideal 
single-channel leads representing a source electrode were 
coupled directly to the styrene molecule at one end of 
the styrene chain, one ideal lead coupling to each of the 



atomic valence p-orbitals of each carbon atom of the ben- 
zene ring of the molecule. 18 ideal single-channel loads 
representing the drain electrode were coupled similarly 
to the styrene molecule at the other end of the molecular 
chain. In this case transport through the styrene chains 
via bands derived from molecular HOMO and LUMO 
orbitals was being investigated so the site energies a, of 
the orbitals making up the ideal leads were chosen lying 
within the molecular HOMO and LUMO energy bands, 
respectively. The nearest neighbor hopping matrix ele- 
ment j3 of the ideal leads was chosen so that the widths 
of the energy bands of the ideal leads were comparable to 
the widths of the energy bands derived from the molec- 
ular HOMO and LUMO energy levels. 



VI. RESULTS: VERTICAL TRANSPORT IN 
THE STM GEOMETRY 

A. Negative sample bias 

Fig. shows representative results for the calculated 
current flowing between the tungsten STM tip and the 
styrene/methyl-styrene molecular chain on silicon that is 
depicted in FigEl The silicon substrate is biased nega- 
tively with respect to the STM tip so that occupied states 
of the molecular chain and substrate are being probed. 
Each curve corresponds to a scan of the STM tip along 
the line of molecules at a constant height above the silicon 
substrate^ The lateral position of the highest carbon 
atom of each molecule is indicated by a dotted vertical 
line labelled m (methyl-styrene) or s (styrene). Curves L 
and L' correspond to a low value of the bias such that the 
Fermi level of the tungsten tip is at -12.0 eV as indicated 
by the arrows labelled L in Fig[S] (b) and (c). Since in 
this case the tungsten Fermi level is well above the upper 
edge of the molecular HOMO band in Fig|S]the current is 
due to electron tunneling through the molecules from the 
silicon substrate states to the STM tip. Curve H is for a 
higher bias with the tip Fermi level at -12.6 eV as indi- 
cated by the arrow labelled H in FigUJb) . In this case the 
tungsten Fermi level is within of the molecular HOMO 
band that is located on the molecular benzene rings and 
mediates transport between the silicon substrate and the 
tip. For curves L and H the tungsten atom at the end of 
the tip is at a height 2 A higher than the highest carbon 
atom of the molecular chain while for curve L' the tip 
atom is 3 A higher than the highest carbon atom. (Note 
that for curve H the current shown has been scaled down 
by a factor of 10 for clarity). 

Curve L is similar to experimental scans along the 
molecular chain at low negative sample bias (FigJ2J (a) 
and (b)), although the experimental scans were taken 
at constant current while the calculation was performed 
at constant tip height^ In curve L each molecule ap- 
pears as a distinct current maximum consistent with 
the experimental data^ and the local current over 
each styrene molecule is somewhat stronger than over 
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a methyl-styrene molecule (as in Fig|5] (a)) despite the 
fact that the methyl-styrene molecule is the taller of the 
two^i Interestingly, as the tip is raised higher above the 
molecular chain at constant (low) bias the current flowing 
through the styrene molecules is predicted to decrease ini- 
tially much faster than that through the methyl-styrene, 
resulting in the significantly weaker styrene current seen 
for curve L'. We attribute this to different behavior of 
weak transport resonances (involving both tip states and 
carbon-silicon interface states) for the styrene molecule 
than for the methyl-styrene molecule that is much closer 
to the tip and thus interacts more strongly with it^ 

As the bias increases, the calculated current through 
the methyl-styrene molecules increases more rapidly than 
through the styrene and the relative amplitude of the 
current corrugation along the molecular chain decreases 
along both the styrene and methyl-styrene parts of the 
chain. Both of these effects (that are clearly visible in 
curve H) are also seen experimentally in Fig|3 (a) and 

(b)ja 

The above theoretical results imply that the observa- 
tion of an experimental STM profile similar to L in Fig. 
El signals that the STM tip Fermi level lies between the 
Si valence band edge at the Si surface and the HOMO 
band of the molecular chain, whereas the observation of 
a profile similar to H signals that the tip Fermi level 
lies below the upper edge of the HOMO band. Since a 
transition from an L type profile to a H type profile is 
in fact observed experimentally as the magnitude of the 
substrate bias is increased (see Fig. EJ), our experimen- 
tal and theoretical findings taken together confirm that 
the molecular HOMO band does in fact lie significantly 
below the silicon valence band edge, as we have already 
suggested in Section llVl based on the nature of the corru- 
gation observed along the molecular chain at low negative 
substrate bias. 

Another interesting aspect of the results shown in Fig. 
Eloccurs at the junction of the methyl-styrene and styrene 
chains: At low bias (curves L and L') the current profile 
of the styrene molecule at the junction is quite similar 
to those of the styrene molecules far from the junction, 
and the same is true of the methyl-styrene molecules. 
However, this is not true at high bias (curve H) where the 
current decreases monotonically from the methyl-styrene 
side of the junction all the way through the location of 
the first molecule on the styrene side; the first styrene 
molecule is not visible as a distinct entity, in marked 
contrast to what is seen near the junction in curve L'. 
We attribute this smearing of the boundary between the 
methyl-styrene and styrene chains in the STM image at 
high bias to derealization of electron states along the 
molecular chain associated with overlapping molecular 
HOMO levels on adjacent molecules. This effect is not 
found in the STM images at low bias (curves L and L') 
where the states outside of the molecular HOMO band 
are responsible for conduction. We will return to this 
topic in our discussion of positive sample bias below and 
its relation to lateral transport along the molecular chain 



will be considered at the end of Section lYlII 

Increased smearing of the boundary between the 
methyl-styrene and styrene chains with increasing bias 
is also seen in the experimental STM data in Fig[3 (a) 
and (b). However, experimentally the apparent height of 
the molecular chain varies on both sides of the junction 
but more so on the methyl-styrene side. The reason for 
this difference between theory and experiment is unclear 
at the present time. 

The calculated current peaks for STM tip scans at low 
bias and constant height in the direction orthogonal to 
that of the scans in Fig. [^exhibit internal structure. This 
can be seen in Fig[7fa) where curve sL (mL) is the result 
of such a scan across the third styrene (methyl-styrene) 
molecule from the junction for the same bias and tip 
height as scan L in FigHjJ With increasing bias voltage 
the internal structure becomes less prominent and the 
scans evolve to a simpler single-peak structure. 



B. Positive sample bias 

Fig. |H| shows the calculated current between the STM 
tip and the molecular chain for some values of positive 
substrate bias. The tip is scanned along the molecular 
chain at constant height^ 2 A higher than the high- 
est carbon atom in each case. The silicon substrate is 
biased positively with respect to the STM tip so that 
unoccupied states of the molecular chain and substrate 
are being probed. Dotted vertical lines again indicate 
positions of the highest carbon atoms of the styrene (s) 
and methyl-styrene molecules. The tungsten Fermi level 
is located below the molecular LUMO band for curves 11 
and 1, slightly above the lower edge of the LUMO band 
for curve e, and well inside the LUMO band for curve 
h, at -9.3 eV, -8.6 eV, -8.39 eV and -8.2 eV respectively, 
as is indicated by arrows in Fig. Efb). At low positive 
bias (curves 11 and 1) electrons tunnel from the tungsten 
tip to the silicon substrate states through the molecules 
outside of the energy range of the molecular LUMO band 
and the current shows a peak over each molecule, as for 
the case of low negative substrate bias (curves L and 
L' of Fig. |HJ). As the bias increases and the tungsten 
Fermi level begins to enter the molecular LUMO band, 
the current develops a minimum over the center of each 
styrene molecule because the molecular LUMO has ir- 
like character with a wave function node in the plane of 
the benzene ringi£& A similar effect also happens over the 
methyl-styrene molecules with increasing bias, but more 
slowly because of the presence of the methyl group be- 
tween the benzene ring and the STM tip: The amplitude 
of the current corrugation over the methyl-styrene chain 
decreases from curve 11 to curve e and then reverses so 
that for curve h the current is depressed over the center 
of each molecule. The amplitude of the current corruga- 
tion also decreases with increasing bias (from curve e to 
curve h) over the styrene chain. At low bias the contrast 
between the methyl-styrene and styrene decreases with 
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decreasing bias (from curve 1 to curve 11) but does not 
reverse as in Fig|SJ 

Interestingly, the current for curve e develops a strong 
overall slope from molecule to molecule extending across 
the entire styrene chain, an effect that is not found for 
curves 11, 1 and h or under negative substrate bias (Fig. 
HjJ). There is also a weaker slope across the methyl- 
styrene chain. This suggests that electronic communi- 
cation between the molecules over a length of several 
inter-molecular spacings is important at this bias volt- 
age and strongly affects the calculated STM image of the 
molecular chain. This will be explored further in Section 

ivm 

FigEJb) shows the calculated tip current for scans 
across the molecular chain at positive substrate bias in 
the direction orthogonal to the scans in Fig. |SJ Curve se 
(me) is the result of such a scan across the third styrene 
(methyl-styrene) molecule from the junction for the same 
bias and tip height as scan e in FiglHl Note the double 
peaked structure in scan se. This double peaked struc- 
ture over styrene evolves into a single peak with increas- 
ing and decreasing bias, while the scan across methyl- 
styrene remains single-peaked. 



VII. LATERAL TRANSPORT AND ITS 
INFLUENCE ON STM IMAGES 

The overlaps between the molecular HOMO and 
LUMO states of adjacent styrene and methyl-styrene 
molecules in the molecular chains on silicon have the 
potential to give rise to bands of electronic states de- 
localized along the molecular chainsAiS However, the 
molecular LUMO (HOMO) energy levels are degenerate 
with states of the silicon conduction (valence) band and 
clearly must hybridize with them to some degree. More- 
over these molecular energy levels correspond to excited 
states of the total system and electrons (holes) occupy- 
ing them must eventually decay to lower energy conduc- 
tion (valence) states of the substrate due to inelastic pro- 
cesses. However, it is still interesting from a theoretical 
perspective to determine how much the hybridization be- 
tween a molecular chain's LUMO and HOMO bands and 
the silicon conduction and valence bands affects electron 
or hole transport along short stretches of the molecu- 
lar chains. For example, if an electron (hole) is injected 
into the LUMO (HOMO) orbital of a molecule belonging 
to the chain and extracted by a probe at another site 
along the molecular chain, is the transport between the 
two sites dominated by conduction through the molecu- 
lar chain or through the substrate or do these two con- 
duction channels make comparable contributions to the 
transport? A related quantity, the characteristic distance 
that an electron or hole injected into a LUMO or HOMO 
state of the molecular chain from an STM tip travels 
along the chain before being absorbed into the silicon 
substrate, may influence the STM image of a molecular 
chain, particularly near a heterojunction. In this Section 



the issue of parallel transport along the molecular chain 
and through silicon substrate over nanoscale distances is 
addressed theoretically, and the results are related qual- 
itatively to the findings regarding the calculated STM 
images of the molecular chains that were presented in 
Section Ell 

The approach used is to calculate the multichannel 
Landauer transmission probabilities for electrons be- 
tween a pair of ideal electrodes coupled to the ends of 
a molecular chain of 8 styrene molecules on a dimer- 
ized hydrogen-terminated silicon substrate. For the cal- 
culations presented here 9 layers of silicon atoms were 
included in the cluster modeling the silicon substrate. 
However, results obtained with a 5 silicon layer model of 
the cluster were similar. The ideal source and drain elec- 
trodes were coupled to the benzene rings of the molecules 
at the two ends of the styrene chain and were represented 
by arrays of ideal single channel leads as described at the 
end of Section |V| The energy bands of the ideal leads 
were taken to be centered at on = —12.7 eV (-8.25 eV) 
for studies of the lateral transport in the energy range 
of the molecular HOMO (LUMO) band; the ideal lead 
hopping parameter (3 was chosen to be -0.25 eV for an 
ideal lead band width of 1 eV. The hamiltonian matrix 
elements coupling the ideal leads to the carbon p-orbitals 
of the molecular benzene rings was chosen to be 0.4/3. 

The calculated Landauer transmission probability 
T(E) (per spin channel) through the system for energies 
corresponding to the LUMO band of the styrene chain is 
shown in Fig. Ela). These results show that the system 
is a fairly good one-dimensional conductor with at least 
two partially-transmitting channels per spin in the en- 
ergy range around the middle of the LUMO band since 
1 < T<2 through much of this energy range. 

To determine whether this strong transmission is oc- 
curring mainly by direct electron transmission from 
molecule to molecule along the chain or mainly through 
the silicon substrate the calculation was repeated with 
all Hamiltonian matrix elements and overlaps between 
atomic orbitals on different molecules switched off so as 
to allow conduction along the chain to proceed only via 
the substrate. (Note that the parameters describing the 
ideal leads and their coupling to the sample were kept 
unchanged.) The result is shown in Fig|5fb). Clearly 
throughout most of the LUMO energy band the transmis- 
sion via the substrate channel (Fig|^b)) is much weaker 
than the direct transmission through the molecular chain 
itself. The transmission through the substrate is strong 
only at a small number of narrow resonances. Further- 
more at the energies at which the strongest of these res- 
onances occur the transmission of the fully coupled sys- 
tem (Fig|5Ja)) shows minima] this means that the reso- 
nant states of the substrate that give rise to the strongest 
transmission peaks in Fig^b) give rise to resonant back- 
scattering (as opposed to contributing to resonant for- 
ward transmission) in the fully coupled system where 
transport is permitted through both the molecular chain 
and substrate. Thus it is quite clear that in the energy 
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range of the LUMO band of the molecular chain, trans- 
port from one end of the molecular chain to the other is 
strongly dominated by direct electron transmission from 
molecule to molecule S 

The corresponding results for transmission from one 
end of the molecular chain to the other in the energy 
range of the molecular HOMO band are shown in FigEJc) 
for the fully coupled system, and in FigHJd) for the sys- 
tem with all Hamiltonian matrix elements and overlaps 
between atomic orbitals on different molecules switched 
off. The molecular chain is still a fairly good one- 
dimensional conductor (T(E) ~ 1 in most of this energy 
range in Fig|5fc)). Transmission through the substrate 
(Fig'Hd)) is again clearly less efficient than through the 
molecular chain although the difference between the two 
is much less stark for the HOMO energy range than for 
the LUMO. 

A simple measure of the overall relative importance 
of direct transmission from molecule to molecule along 
the molecular chain relative to that via the parallel path 
through the substrate is the ratio R of the average trans- 
mission for the fully coupled system to the average trans- 
mission with the molecule-to-molecule coupling switched 
off. For the 8 styrene molecule chain R ~ 30 for the 
molecular LUMO band while for the molecular HOMO 
band R ~ 4. 

For ideal source and drain leads coupled to opposite 
ends of the substrate instead of to the molecular benzene 
rings at the ends of the styrene chain, the calculated aver- 
age transmission probabilities in the HOMO and LUMO 
energy ranges were found to be similar in magnitude. 
Thus the above difference in R between the HOMO and 
LUMO bands of the molecular chain is mainly due to 
weaker hybridization of the substrate orbitals with the 
molecular LUMO states and than with the molecular 
HOMO. 

Since the above calculations indicate that states of the 
LUMO band of the molecular chain hybridize much less 
with the states of the silicon substrate than the states of 
the molecular HOMO band do, it is reasonable to expect 
it to be easier for STM measurements to detect a distinc- 
tive signature of electronic states confined to the molecu- 
lar chain and extending over a significant distance along 
the molecular chain in the LUMO energy range than in 
the HOMO range. The results of the calculations of ver- 
tical transport in the STM geometry presented in Section 
IVII are consistent with this: In particular, one of the the- 
oretical STM current plots (curve e) in Fig. [S] shows 
an overall decline in the STM current as one moves away 
from the center of the molecular chain in either direction; 
the decline extends all the way to the end of the molec- 
ular chain in either direction but is especially strong on 
the styrene side of the chain. At the bias voltage cor- 
responding to this current plot, the Fermi level of the 
STM tip is slightly above the lower edge of the LUMO 
band of the molecular chain and a single transmission 
resonance due to lowest electronic eigenstate belonging 
to the LUMO band makes the dominant contribution to 



the current. The electron probability distribution in this 
eigenstate is peaked at the center of the molecular chain 
and declines strongly towards both ends of the molecu- 
lar chain. This explains qualitatively why the calculated 
STM current (curve e, Fig. |SJ is stronger at the center 
of the chain than at its ends. (However, to predict the 
current profile along the chain quantitatively knowledge 
of the LUMO eigenfunction alone is not sufficient and 
a complete transport calculation such as is described in 
Section[V]is necessary.) By contrast, the calculated STM 
current for the case where the HOMO band of the molec- 
ular chain is being probed does not show such pronounced 
trends across more than a single molecule, even when the 
tungsten Fermi level is close to the edge of the HOMO 
band, just as one might expect qualitatively for a molec- 
ular band that mixes more freely with states belonging 
to the substrate. 



VIII. END EFFECTS 

The model geometries of the molecular chains studied 
theoretically in the preceding sections were constructed 
(see Section llllfl to be representative of regions of the 
chain far from its ends. This was done by constraining 
the average spacing between the tops of the molecules 
to be similar to that observed experimentally near the 
center of the molecular chain. If this constraint is not 
imposed, molecules at the ends of the chain tilt more 
strongly towards the silicon substrate. This is seen ex- 
perimentally in Fig. |2t a) and (b) , where the apparent dif- 
ference in height between the end and central molecules 
of the chain at low negative sample bias is more than 
an Angstrom. However, the observed height contrast be- 
tween the ends and center of the chain in Fig. Ufa) and 
(b) decreases markedly as the magnitude of the bias in- 
creases. Thus considerations other than the geometrical 
height profile of the molecular chain must also be taken 
into account to explain this data. 

The gaps between the molecular HOMO and LUMO 
bands of stacks of parallel-oriented aromatic molecules 
should increase with increasing separation between 
the molecules because the coupling between adjacent 
molecules decreases as their separation increases^ This 
suggests the following scenario as a simple explanation 
of the bias-dependence of the STM profiles of the molec- 
ular chains that we have observed: The tilting apart of 
the molecules near the ends of the chains causes the lo- 
cal HOMO-LUMO band gap of the molecular chain to 
increase towards the ends of the chain. In constant cur- 
rent STM images at low bias where the Fermi level of 
the tip is in the molecular HOMO-LUMO band gap, this 
increased band gap causes the molecular chain to appear 
lower at its ends relative to its center than the geomet- 
rical difference in height alone would require. At larger 
bias, when the tip Fermi level has fallen below the (local) 
top of the molecular HOMO band throughout the molec- 
ular chain and all of the molecules have "turned on," 
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the end molecules appear less diminished. It appears 
therefore that the variable HOMO-LUMO gap scenario 
could account qualitatively for the position and voltage 
dependent height variations observed. However, a more 
thorough analysis shows this model to be incomplete; it 
does not take into account three aspects of this system 
that we show below to be crucial. They are: (1) The 
molecules bond to equally spaced Si dimers on the (001) 
silicon surface. Thus although the spacing of the tops of 
the molecules in STM images varies along the chain due 
to the varying tilt of the molecules from the surface nor- 
mal, the spacing between the bottoms of the molecules 
is fixed and thus, as will be seen below, even the most 
strongly tilted molecule at the end of the chain couples 
strongly to the rest of the chain. (2) The coupling of 
molecular HOMO and LUMO orbitals to the substrate 
(to which vertical transport is very sensitive) increases 
strongly along the chain as its end is approached be- 
cause of the increasing tilt of the molecules towards the 
substrate^ (3) At low bias where the STM tip Fermi 
level is in the molecular HOMO-LUMO gap the trans- 
mission of electrons between the substrate and tip is me- 
diated by resonant states at the C-Si interface rather than 
the molecular HOMO and LUMO orbitals, as is discussed 
in Section IVT1 i.e., a different transport mechanism that 
is not included in the pure variable band gap scenario 
outlined above operates in the low bias regime. 

For simplicity, we restrict our attention to styrene 
chains. Our results for the (vertical) STM current over 
a chain of 8 styrene molecules on a hydrogen-terminated 
dimerized Si(100) substrate for negative sample bias are 
shown in Fig. ^3 The silicon substrate, tungsten STM 
tip and ideal source and drain leads are as in Section lv"fl 
However, the structure of the molecular chain is different: 
Here the final AMI relaxation of the chain differed from 
that described in Section lTTll in that the styrene molecule 
at the left end of the chain was constrained to stand 
upright on the substrate (with a small tilt to the right) 
but the other molecules were not constrained in this way. 
Thus the structure of this model styrene chain near its 
left end resembles that in the interior of long molecular 
chain (as for the model structures considered in Section 
IVI|) while its right end should be similar to that of the 
real end of a molecular chain. The currents shown in 
Fig. ^]]are calculated for STM tip trajectories that pass 
directly over the 4th carbon atom of the benzene ring of 
each styrene molecule, i.e., the carbon atom furthest from 
the carbon chain that anchors the benzene ring to the sil- 
icon. [The locations of these carbon atoms are shown by 
dotted vertical lines in Fig. E3] This 4th carbon atom 
for the molecule at the right end of the model styrene 
chain is located 1.12 Acloser to the Si substrate than the 
highest carbon atom of the chain^Siffi a height difference 
similar to the apparent height contrast between the end 
and interior molecules of the chain observed experimen- 
tally at low bias; see Fig[2Ja,b). The current plots L (H) 
in Fig. EH are for the same low (high) values of the neg- 
ative sample bias as the plots L (H) in Fig. for L (H) 



the STM tip Fermi level is above (within) the molecular 
HOMO band as in Fig. |3J However, in Fig. ^|the cur- 
rents labelled H have been scaled up by a factor of 5 for 
clarity while that in Fig. [5] was scaled down by a factor 
of 10. The solid curves L and H in Fig. ^] show the 
calculated STM current for scans at constant tip height 
above the Si substrate (the same height as for curves L 
and H in Fig. [SJ. The dashed curves, however, show 
calculated currents for tip trajectories which follow the 
sloping height profile of the relaxed molecular chain (thus 
at varying height from the underlying silicon surface). 61 

For the scans at constant height above the silicon sur- 
face (the solid curves) the calculated current over the 
tilted molecule at the right end of the chain is lower than 
that near the left end, by a factor of more than 30 for the 
low bias case L and a factor of ~ 10 for the higher bias 
case H. These results agree qualitatively with our exper- 
imental finding that the height of the molecules appears 
lower at the end of the chain and that the magnitude of 
the lowering decreases as the magnitude of the bias in- 
creases. However for the lower bias scan at a constant 
vertical distance above the tops of the molecules (dashed 
curve L) the calculated current does not decrease near 
the right end of the molecular chain; in fact it is slightly 
higher over the molecule that ends the chain than over the 
rest of the chain. By contrast, for the simple position- 
dependent band gap model outlined above to apply in 
isolation the apparent height of the molecules at the end 
of the chain must be substantially lower than their ge- 
ometrical height. Thus our calculations do not support 
a purely position-dependent band gap scenario. They 
suggest instead that at low negative substrate bias in 
the constant current STM mode the experimentally mea- 
sured differences in the height of the tip from styrene 
molecule to styrene molecule along the chain should, to a 
first approximation, be equal to the differences in the ge- 
ometrical heights of the tops of the respective molecules. 
This interpretation is further supported by the approxi- 
mate agreement between the apparent height difference 
in experimental data in Fig. |2Ja,b) between the styrene 
molecules near the center and end of the chain at low bias 
and the result of our calculation of the structure at the 
end of the molecular chain: In each case the top of the 
molecule at the end of the chain is lower than the tops of 
the molecules far from the ends by somewhat more than 

l A. 

We explain the bias dependence of the apparent height 
of the molecules at the end of the chain relative to those 
far from the end as follows: Because the molecules near 
the end of the chain tilt more strongly towards the sub- 
strate than the other molecules of the chain do, the 
molecular HOMO orbitals of the end molecules (being 
mainly on the benzene rings) couple much more strongly 
to the substrate than the HOMO orbitals of the other 
molecules do. This stronger coupling leads to stronger 
electron transmission from the substrate to the STM tip 
through the end molecules of the chain than through 
the other molecules at the higher values of the negative 
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substrate bias for which the Fermi level of the tung- 
sten tip falls below the top of the molecular HOMO 
band so that transport via the HOMO orbitals domi- 
nates the conduction^ Thus for the dashed curve H in 
Fig. where the tip scans along the chain at same 
vertical distance above each molecule, the current grows 
quite significantly as the end of the chain is approached. 
Consequently the decrease in the apparent height of the 
molecules as the end of the chain is approached in con- 
stant current STM mode is not as strong as it is at lower 
bias where the molecular HOMO states do not mediate 
the transport^ 

This mechanism is clarified further by the inset of Fig. 
ITU1 Here the Landauer transmission T(E, V) between the 
substrate and STM tip is plotted as a function of elec- 
tron energy E for the same (higher) value of the applied 
bias V as for curves H of Fig. ^3 For each transmis- 
sion curve plotted the STM tip is located 3.497 Aabove 
the 4th carbon atom of the benzene ring of the nearest 
molecule (as for the dashed curve H of Fig. llOjl . For 
reference the arrow indicates the HOMO energy of an 
ethyl-benzene molecule in free space. The solid line is the 
calculated transmission with the tip above the molecule 
at the right end of the molecular chain (with the other 
molecules of the chain also present on the substrate). The 
dashed line shows the transmission with the tip above the 
same molecule (with the same tilted geometry) but with 
the other molecules of the chain removed from the sys- 
tem and replaced by hydrogen atoms. The dot-dashed 
curve is the average transmission for tip positions over 
the 5 molecules closest to the left end of the molecular 
chain. The nearly equally spaced peaks of the dot-dashed 
curve arise from states of the molecular HOMO band that 
is formed due to the coupling between the HOMO or- 
bitals of the molecules making up the chain** The much 
greater strength of the transmission for the tip above the 
molecule at the right end of the chain (solid curve) is 
due to the stronger coupling of that molecule's HOMO 
orbital to the substrate, as discussed above. It results in 
the stronger current when the tip is over the right end 
molecule than over the rest of the chain, as is seen in 
the dashed curve H. The stronger coupling of the right 
end molecule's HOMO to the substrate also results in 
the strong broadening of the the associated transmission 
peaks (solid and dashed lines of the inset); the HOMO 
transmission peaks (not shown) for isolated molecules on 
the substrate with the upright geometries of molecules 
far from the end of the chain are much narrower, with 
widths similar to those of the individual peaks of the 
dot-dashed curve of the inset. 

Although, as we have shown, the enhanced coupling 
of the molecular HOMO at the end of the chain to the 
substrate plays a crucial role in the bias dependence of 
the STM profile of the end of the chain, the coupling 
of the end molecule to the rest of the molecular chain 
remains strong despite the tilting of the end molecule 
towards the substrate. This is evident from a comparison 
of the transmission peaks for the end molecule with the 



rest of the molecular chain present and absent from the 
system: The coupling of the end molecule to the rest of 
the chain results in the large splitting of the transmission 
peak that occurs when the rest of the molecular chain is 
present (solid curve of the inset) and does not occur when 
the rest of the chain is absent (dashed curve) 

Our results for positive substrate bias are shown in Fig. 
1111 The chain geometry and the meaning of the solid and 
dashed plots are as in the negative bias case (Fig. ITU)l . 
As in Fig. [HI for plots 1, e and h the STM tip Fermi level 
is below, slightly above and well above the lower edge of 
the molecular chain's LUMO band (see Fig. |SJ), but in 
Fig. ^2 the currents 1 have been scaled down by a factor 
of 10. The molecules at the end of the chain are again 
predicted to appear lower in STM images than those near 
the center, and this lowering becomes less pronounced as 
the magnitude of the applied bias voltage increases^ 
The mechanism of the bias dependence is similar to that 
for negative substrate bias: Because of the tilting of the 
molecules at the end of the chain, their LUMO orbitals 
couple more strongly to the substrate, and this results 
in stronger and broader transmission peaks (and an en- 
hanced STM current) over the end molecules when the 
applied bias is large enough for the LUMO orbitals to 
participate in conduction. However, for plots e and h of 
Fig. II II the local current maxima occur over or near the 
tops of the molecules near the end of the chain while far 
from the end current minima coincide with the molecular 
positions. Because of the long distance that the influence 
of the physical (right) end of the molecular chain extends 
along the chain in plot e, experimentally distinguishing 
effects related to the presence of a hetero junction from 
such end effects in empty-state STM imaging will require 
long molecular chains and careful analysis. 



IX. CONCLUSIONS 

We have synthesized heterojunctions of lines of styrene 
and methyl-styrene molecules on H-terminated (100) sili- 
con substrates and have studied these systems with scan- 
ning tunneling microscopy as well as theoretically. We 
have developed a tight binding theory based on the ex- 
tended Hiickel model (modified so as to provide a good 
description of the electronic structure of tungsten and 
silicon crystals as well as of the molecules) that together 
with Green's function and Landauer transport techniques 
is able to account for key features of the experimental 
filled-state STM images of these molecular systems (in- 
cluding both the heterojunctions and the regions near 
the ends of the molecular chains) and of their variation 
with applied bias voltage. This model has also been used 
to explore theoretically the lateral transport properties 
of short chains of styrene molecules on silicon and it has 
been found that electron transmission along such chains 
should occur primarily from molecule to molecule rather 
than through the parallel substrate channel, especially 
in the energy range of the molecular LUMO band. The 
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calculations also suggest that phenomena due to elec- 
tron wave functions extending from molecule to molecule 
(such as states localized around the heterojunction of a 
pair of molecular wires) may be observable experimen- 
tally by STM imaging. 
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FIG. 1: (Color on line). Constant current filled-state STM imaging of a methyl-styrene/styrene heterostructure on the H:Si(100) 
surface, a, Growth of the first line segment of methyl-styrene (8L exposure). The chemically reactive silicon radical which 
initiated the line growth (not shown) is now located at the end of the CH3-styrene segment (indicated by white arrow). This 
silicon radical will serve as the initiation site for the growth of the second line segment, b, Following a 78 L exposure of styrene, 
the methyl-styrene segment in panel (a) has been extended to form a molecular heterostructure. The white wedge indicates 
the location of the heterointerface. At V s = -3.0 V, the methyl-styrene line segment images higher than the styrene segment, 
c, V s = -2.4 V, the methyl-styrene and styrene line segments image with similar height, d, V s = -1.8 V, the methyl-styrene 
and styrene line segments reveal different molecular contrast. Inset: Blue and red arrows show locations of topographic cross- 
sections presented in Figures|!2](a) and (b), respectively. Images were acquired using a constant tunnel current of 60 pA. Image 
areas: 8.5 nmx8.5 nm. 

FIG. 2: (Color on line). Topographic cross-sections extracted from STM image data of the methyl-styrene/ styrene het- 
erostructure presented in Figure 1. a, Topographic cross-sections taken along the chemical attachment points between the 
heterostructure and the underlying dimer row - along direction of blue arrow in Figure 0(d). b, Topographic cross-sections 
taken to the right of the trench separating the underlying dimer rows - along direction of red arrow in Figure (d) . Cross- 
sections for a range of sample biases (-3.0 V to -1.8 V) are shown in each panel. Black arrows in panels (a) and (b) indicate 
the position of the interfacial methyl-styrene molecule. 

FIG. 3: (Color on line). Model structure of a heterojunction between methyl-styrene and styrene chains (each of four molecules) 
on a hydrogen-passivated silicon cluster representing a dimerized H-terminated Si (001) surface. Carbon atoms are black, silicon 
is blue, hydrogen is white. Inset: Another view of the molecules and nearby substrate atoms. 

FIG. 4: (Color on line). Geometry of a methyl-styrene molecule and styrene molecule relaxed using density functional theory 
on a small dimerized Si (001) cluster. Notation as in Fig[3] For significance of the numbering see text. 

FIG. 5: Histogram of the silicon (a), styrene carbon (b), methyl-styrene carbon (c) and hydrogen (d) content of the electronic 
eigenstates of the eight-molecule styrene/methyl-styrene chain on the hydrogen passivated silicon cluster depicted in Fig|3] The 
arrows labelled H and L indicate the positions of the tungsten Fermi level relative to the molecular energy spectrum at the 
negative substrate biases corresponding to current curves H and L,L' in FigHJ] Arrows 11, 1, e and h show the tungsten Fermi 
level at the positive substrate biases corresponding to current curves 11, 1, e and h in Fig|5] 

FIG. 6: Calculated current I (arbitrary units) flowing between the tungsten STM tip and the styrene/methyl-styrene molecular 
chain on silicon at some negative substrate biases vs. STM tip position along the chain. The current for curve H has been 
scaled down by a factor 10. Dotted vertical lines labelled m (s) indicate the locations of the highest carbon atoms of the 
methyl-styrene (styrene) molecules. The scans are at constant tip height, 2 A (3 A) higher than the highest carbon atom of 
the chain for curves H and L (curve L'). L and L' are at a lower bias, H at a higher bias. The tungsten Fermi level at each 
bias is indicated by an arrow in Fig|S] see text. 

FIG. 7: Calculated current between the tungsten STM tip and the styrene/methyl-styrene molecular chain on silicon vs. STM 
tip position scanned across the chain, in the direction perpendicular to that of the scans in Fig|S|and|H| Notation as in Figlol 
The scans are at constant tip height of 2 A higher than the highest carbon atom of the molecular chain, (a) Scan sL (mL) is 
over the third styrene (methyl-styrene) molecule from the junction for the same bias as scan L in Fig|S| (b) Scan se (me) is 
over the third styrene (methyl-styrene) molecule from the junction for the same bias as scan e in Fig|S] Each scan passes over 
the highest carbon atom of the respective styrene or methyl-styrene molecule and crosses the locus of the corresponding scan 
of Fig|S|or|H|at position = 0. 

FIG. 8: Calculated current between the tungsten STM tip and the styrene/methyl-styrene molecular chain on silicon at some 
positive substrate biases vs. STM tip position along the chain. Notation as in FigHJbut the scale is the same for all curves. 
The scans are at constant tip height of 2 A higher than the highest carbon atom of the molecular chain. The magnitude of the 
bias voltage increases from curve 11 to 1 to e to h (see text). The tungsten Fermi level at each bias is indicated by an arrow in 
FigEl 

FIG. 9: Calculated Landauer transmission probability T(E) (per spin channel) through a molecular chain of 8 styrene molecules 
on a hydrogen-terminated silicon cluster in the energy range of the molecular LUMO band (a) and of the molecular HOMO 
band (c). Transmission in the same energy ranges with all Hamiltonian matrix elements and overlaps between orbitals of 
different molecules switched off (so that transport between the terminal molecules must proceed via the substrate) is shown in 
plots (b) and (d), respectively. 
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FIG. 10: Calculated current I (arbitrary units) flowing between the tungsten STM tip and a styrene molecular chain on silicon 
at some negative substrate biases vs. STM tip position along the chain projected onto an axis parallel to the Si dimer rows 
of the substrate. The tungsten Fermi level at each bias (L or H) is indicated by an arrow in FigQ>] The styrene molecule 
at the left end of the chain has been constrained to stand approximately upright in a geometry representative of the region 
near the center of a molecular chain while the other molecules of the chain have been allowed to relax so that the molecular 
geometries at right end of the chain are representative of the physical end of a molecular chain on silicon. (See text). The 
currents for curves H has been scaled up by a factor 5. Dotted vertical lines labelled s indicate the lateral locations of the 
carbon atoms of the styrene molecules that are furthest from the carbon chains anchoring the molecules to the silicon. The 
dashed curves are for scans in which the tip passes over each styrene molecule at an equal vertical distance (3.497 A) above 
this carbon atom. The solid curves are for scans in which the tip is kept at a constant height above the silicon substrate. Inset: 
The solid curve shows the Landauer transmission probability between tip and substrate vs. electron energy in eV with the 
STM tip 3.497 Aabove the styrene molecule at the right end of the chain for the same applied bias as curves H. Dashed curve: 
Same except that all molecules other than the molecule at the right end of the chain are replaced with H atoms passivating the 
Si substrate. Dot-dashed curve: Landauer transmission probability averaged over tip positions 3.497 Aabove the five styrene 
molecules nearest the left end of the molecular chain. 



FIG. 11: Calculated current I flowing between the tungsten STM tip and a styrene molecular chain on silicon at some positive 
substrate biases vs. STM tip position. The tungsten Fermi level at each bias (1, e and h) is indicated by an arrow in Figl^l 
Notation as in Fig. ^| 
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Fig. 6, Kirczenow, Piva, and Wolkow. 
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Fig. 8, Kirczenow, Piva, and Wolkow. 
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Fig. 1 1 , Kirczenow, Piva, and Wolkow. 



